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Abstract 

We study the detectability of primordial metal-free stars, both through direct 
searches for their emission, as well as searches for the resulting supernovae. We 
show that enrichment is a local process that takes place over an extended redshift 
range. While the duration of the transition from a metal-free to an enriched universe 
depends on several unknown factors, in all models late-forming metal-free stars are 
found in 10^'^ — 10^ '^ Mq objects, which are just large enough to cool within a 
Hubble time, but small enough to not be clustered near areas of previous star for- 
mation. We discuss the observational properties of these objects, some of which may 
have already been detected in ongoing surveys of high-redshift Lyman-a emitters. 

If metal- free stars have masses 140 ^ ^ 260 Mq, they are expected to end 
their lives as pair-production supernovae (PPSNe), in which an electron-positron 
pair-production instability triggers explosive nuclear burning. Using the implicit 
hydrodynamics code KEPLER, we calculate a set of PPSNe light curves that allows 
us to assess observational strategies for finding these objects. The peak luminosities 
of typical PPSNe are only slightly greater than those of Type la supernovae, but they 
remain bright much longer (~ 1 year) and have hydrogen lines. Ongoing supernova 
searches may soon be able to place stringent limits on the fraction of very massive 
stars out to ^ ~ 2. The planned Joint Dark Energy Mission satellite will be able to 
extend these limits out to z ~ 6. 
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1 Introduction 



A long time ago tliere were no metals. The primordial fireball generated a 
whirlwind of neutrinos, an ocean of thermal photons, and a mountain of light 
nuclei, but failed to produce elements heavier than lithium. In fact, it was not 
until several hundred million years later that the universe was first enriched 
with heavy elements. 

Yet, today stellar nucleosynthetic products are everywhere. The lowest metal- 
licity galaxies known are enriched to substantial values of ~ 0.02^© (0); the 
lowest-density regions of the intergalactic medium (IGM) appear to be en- 
riched out to the highest redshifts probed and although direct obser- 
vations of metal-poor halos stars have uncovered a number unusual abundance 
patterns not a single metal-free star has ever been detected. 



Thus a clear signature of the transition from a metal-free to an enriched uni- 
verse has so-far remained unobserved. Furthermore, there are good reasons to 
believe that this transition may have marked a drastic change in the character- 
istics of stars. Recent theoretical studies have shown that collapsiiig primordial 
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clouds fragmented into clumps with typical masses ~ 10 M, 
that accretion onto proto-stellar cores within these clouds was very efficient 
(3 3)- The implication is that the initial mass function (IMF) of primordial 



stars may have been biased to higher masses than observed today. 



2 The Epoch of Metal- Free Stars 



The work discussed in this section and in S3 is described in further detail 



in As the detailed shape of the metal-free IMF is highly uncertain, we 
adopted a wide range of theoretical models that are based on studies of pri- 
mordial gas cooling and fragmentation BBS IS)- In these investigations, 
the minimum fragment mass was comparable to the Jeans mass at the tem- 
perature and density at which molecular hydrogen levels start to be populated 
according to LTE, which is ~ 10^ M©. 

With these studies in mind, we parameterized the IMF with a range of mod- 
els in which the IMF is a Gaussian of the form -7= — ^-iM-Mc) /so-^ where 

V ZncTc 



100 Mq < Mc < 1000 Mq and < ac < a^?^^. Finally, we also explored 
a "null hypothesis" in which the distribution metal-free stars is given by a 
Salpeter IMF. These models can be parameterized in turn by a single quan- 
tity, E^^^^ the ' 'energy input per unit gas mass," which is the product of the 
fraction of gas mass that goes into stars (/^), the number of supernovae per 
M0 of stars (Ngjj), the kinetic energy input per SN in units of 10^^ ergs {E^^/), 
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and the fraction of the SN kinetic energy that gets channeled into an outflow, 
dispersing the metals (/„) outside the host galaxy. Thus E^^^ = f^NgjjE^^^f^ 
has units of 10^^ ergs Mq"^. 



Zero-metallicity stars with initial masses 140 Mq ^ ^ 260 Mq are expected 
to end their lives as pair-production supernovae PPSNe, in which an electron- 
positron pair-production instability triggers explosive nuclear burning (14; ish . 
In this case, the kinetic energy per SN is ~ 50 x 10^^ ergs, and the ejected 
metal mass is ~ 100 (jish. In conventional core-collapse supernovae, the 
energy per SN is about 10^^ ergs and roughly 2 Mq of metal are ejected. Thus 
regardless of the IMF, the total mass in metals is roughly 2 Mq per 10^^ ergs 
of kinetic energy input. 



In other words, despite the many uncertainties surrounding metal-free stars, 
the total energy in each outflow and the total metal mass dispersed are closely 
related. Using a simple 1-D outflow model (0) to study the metallicity evolu- 
tion within each dispersal event, we have found that for the full range of pa- 
rameters considered, the mean metallicity of outflows from primordial objects 
is well above the critical threshold {Zcr ~ 10~'^Zq) that marks the formation 
of normal stars ([lol). This means that primordial star formation continued well 
past the time at which the average IGM metallicity reached Z^r, as large areas 
of the IGM remained pristine while others were enriched to many times this 
value. Enrichment was a local process, and the history of metal-free formation 
is likely to have been determined almost exclusively by the efficiency of metal 
ejection, parameterized by -Ef^^- 



To track the dispersal of the first ejected metals, we made use of the analytical 
formalism described in ([iTi l [see also ([lih] coupled with our 1-D outfiow model. 
The resulting evolution of the star formation rate (SFR) density and the mean 
cosmic metallicity are shown in Figure 1. Here E^^^ ranges from the 10~^ value 
appropriate for a Salpeter IMF to the extreme 10~^'^ value appropriate for an 
IMF composed almost solely of ~ 200 Mq stars. 



In all models, while the peak of primordial star formation occurs at z ^ 10, 
such stars continue to contribute appreciably to the SFR density at much lower 
redshifts. This is true even though the mean IGM metallicity has moved well 
past the critical transition metallicity. Furthermore, for many models, sub- 
stantial primordial star formation continues well into the observable redshift 
range, contributing to as much as 10% of the SFR at ^ = 5. In all cases these 
objects tend to be in the lO'''^ Mq — 10^'° Mq mass range, just large enough to 
cool within a Hubble time, but small enough that they are not clustered near 
areas of previous star formation. 
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Fig. 1. Star formation rate densities for primordial and PopII/I objects. All models 
assume f:^^^ = 0.1, = 0.1, and = 0.1. Left : SFR density per cubic comoving 
Mpc for primordial objects (dashed lines) and PopII/I objects (solid lines). We 
consider models in which, from top to bottom, log]^o(-^g^^) = —4.0 —3.5, —3.0, and 
—2.5. Right: IGM average metallicity resulting from primordial (dashed lines) and 
PopII/I (solid lines) objects. Again, log^Q(£'g^^) ranges from —4.0 to —2.5 from top 
to bottom. 

3 Direct Detections of Primordial Objects 



Employing a simple model of burst-mode star formation, we have developed 
a rough characterization of these directly observable primordial objects. As 
metal-free stars are powerful Lya line emitters lifl ). it is natural to use 
this indicator as a first step in any search for primordial objects. The Lyct lu- 
minosity, Lq, is directly related to the corresponding hydrogen ionizing photon 
rate, Q{H), by L„ = cl(1 - fesc)Q{H), where cl = 1.04 x 10"^^ erg and fesc 
is the escape fraction of ionizing photons from the galaxy. For this uncertain 
factor we adopt an educated guess of fesc = 0.2, which is based on a compi- 



lation of theoretical and observational results ()2l[ ). We c omp ute Q{H) from 
evolutionary stellar models: either from STARBURST99 ()22n in the Salpeter 
IMF case, or from the models in (jiol ) in the case of very massive stars. 

In Figure 2 the isocontours in the Lya luminosity-redshift plane indicate the 
probability to find a galaxy hosting metal-free stars in a given sample of emit- 
ters, again for various values of E^^^ . In the left panel of Figure 2, we consider 
a model in which metal-free stars form according to a Salpeter IMF. Although 
in this case the Lya emission from these stars is slightly brighter than from 
their Popll/1 counterparts (fl^, this difference is not dramatic. However, as 
outflows are relatively weak {E^^^ ~ 10~^), and thus ~ 10% of 2; = 5 Lya 
emitters with observed fluxes above 1.5 x 10"^^ ergs cm~^ s~^ are likely be 
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Fig. 2. Fraction of primordial objects as a function of Lya luminosity and redshift. 
Isocontours of fractions > 10~^, 10~^'^, 10~^, and lO"'''^ are shown. The left panels 
correspond to a model with a Salpeter IMF, while the right panels correspond to a 
model in which primordial stars are ~ 200 Mq. Each panel is labeled by its assumed 
E^^^ value. For reference, the dashed line gives the luminosity corresponding to an 
observed flux of 1.5 x 10~^^ ergs cm~^ s~^, and the points correspond to a number 
of recent measurements [See ((ill) for details]. 

made up of metal-free stars. 

In the right panel of Figure 2, we consider a model in which the metal-free 
IMF is peaked at ~ 200 Mq. In this case the short but bright evolution of such 
stars boosts their luminosities by over an order of magnitude. However, the 
ejection of the resulting metals is dramatically more efficient [E^^^ ~ 10~^), 
and thus again ~ 10% of 2; = 5 Lya emitters are metal-free. 

Thus it is likely that a number of metal-free objects are lurking at the limit 
of present surveys. The goal then is to find a clear-cut feature that uniquely 
identifies a given object as metal-free. The most simple such criterion is the 
equivalent width of Lya, which should be higher in primordial stars than 
in a normal population, as such stars burn hotter and generate more ion- 
izing photons (jioi ). particularly in the very- massive case ^xh . In fact, such 
large equivalent widths have indeed been seen in the observed populations of 
high-redshifts Lya emitters (jiih. although this can also be the result of other 



mechanisms, such as scattering in a clumpy, dusty medium (|24l). 



A much less ambiguous indicator has been identified by ()25[). who noted that 
primordial stars produce large He HI regions, which may emit detectable He 
II recombination emission. The recombination lines at A1640 A(n = 3 — > 2), 
A3203 A(n = 5^ 3), and A4686 A(n = 4 — >■ 3) are particularly attractive 
for this purpose because they suffer minimal effects of scattering by gas and 
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decreasing attenuation by intervening dust. For star formation rates of 20 
Mq yr~^ and 5 Mq yr~^ in a given galaxy, the A1640 flux is detectable out to 
z ~ 5 at the sensitivity level of current surveys, although the flux for A4686 is 
7.1 times lower. While these lines also arise from Wolf-Rayet stars, they would 
be accompanied by strong metal lines in this case (27). Metal-free stars, on the 
other hand, are uniquely intense sources of pure He-II emission. Note that this 
spectroscopic test is particularly difficult, however, as metal-free stars should 
only be found in a small fraction of ^ = 5 objects, and composite spectra 
are likely to wipe out such a signal. Thus the observed absence of the 1640 
doublet in the combined spectra from the Large Area Lyman Alpha survey 
(jH) is consistent with all but the most optimistic of the models in Fig 2. 



4 Pair-Production Supernovae 



The work discussed in this section is described in further detail in (|26l). If 
primordial stars were indeed very massive, the resulting pair-production su- 
pernovae may be directly detectable. This approach has the advantage of con- 
straining the total rate density of very massive star formation as a function of 
redshift, regardless of how they are distributed in galaxies of different sizes. 

Previously, only very limited or approximate models of PPSNe have been 
available in th^ literature 0HE1E3. Here we use the KEPLER code 
(jisl ) to directly construct a suite of PPSN lightcurves that addresses the range 
of theoretical possibilities and accounts for the detailed physical processes 
involved up until the point at which the SN becomes optically thin. These 
results are shown in Figure 3. 

For most of their lifetimes, the effective temperatures of PPSNe are just above 
the ~ lO'^'^K recombination temperature of hydrogen, which corresponds to a 
peak black-body wavelength f« 8000 A. The most important factors in deter- 
mining the features in the luminosity evolution of these SNe are the mass of 
the progenitor star and the efficiency of dredge-up of carbon from the core into 
the envelope. In general, increasing the mass leads to greater ^^Ni production, 
which boosts the late time SN luminosity, through radioactive decay. Mixing, 
on the other hand, has two major effects: it increases the opacity in the en- 
velope, leading to a red giant phase that increases the early-time luminosity; 
and it decreases the mass of the He core, consequently decreasing the ^^Ni 
mass and lowering the late-time luminosity. 

Despite these uncertain factors, PPSNe in general can be characterized as SNe 
that have all three key features: (1) peak magnitudes that are brighter than 
Type II SNe and comparable or slightly brighter than typical SNe Type la; 
(2) very long decay times ~ 1 year, which result from the large initial radii 
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Fig. 3. Comparison of light curves of a SN Type la, a SN Type ITP, a bright SN Type 
II-L, and PPSNe models with varying progenitor masses. In all cases the solid lines 
are absolute V-band AB magnitudes, the dot-dashed lines are the absolute B-band 
AB magnitudes, and the dashed lines are the absolute U-band AB magnitudes. The 
model names refer to the mass of the progenitor star (in units of Mq) and the weak 
(W), intermediate (I), or strong (S) level of mixing by convective overshoot. 

and large masses of material involved in the explosion; and (3) the presence of 
hydrogen lines, which are caused by the outer envelope. Note, however, that 
only this last feature is present in a// cases, and in fact, the lowest mass PPSN 
models we constructed have lightcurves that are remarkably similar to those 
of SN Type II. 



Accounting for redshifting, time dilation, and appropriate luminosity factors, 
we used these lightcurves to relate the overall very massive SFR density to the 
number of PPSNe detectable in current and planned supernova searches. As 
the metal-free SFR density is unknown, we adopted two simple models, which 
can be easily shifted for comparison with literature models, including our 
models in Fig. 1. In the first model, we assume that metal-free star formation 
occurs at a constant rate density, which we take to be 10~^ M0yr~^Mpc~^, 
independent of redshift. In the second case, we assume that at all redshifts 
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Fig. 4. Number of PPSNe per square degree per unit redshift above a fixed 
broad-band magnitude, assuming 0.001 M0 yr-i Mpc'^ (solid lines) or 1% of the 
observed SFR density (dashed lines). Top: The Jab = 27.7 limit taken in these pan- 
els corresponds to a single scan in the planned deep-field SNAP survey, which would 
cover 7.5 deg^ (indicated by the upper dotted lines). A similar magnitude limit with 
the same coverage would be obtained with 2 months of grism data from the Destiny 
mission. The full planned 700 deg^ wide-field SNAP survey (lower dotted lines) also 
has the same limiting magnitude. Center: The Jab = 30.3 limit taken in these pan- 
els is that of the deep-field SNAP survey, which would be able to place constraints 
on PPSN out to 2 > 5. Bottom: Curves corresponding to the JEDI mission, with 
a magnitude limit of /Cab = 28.7 and a coverage of over 24 deg^ with 2 months of 
data. At the highest redshifts, this search does substantially better than the fainter, 
but bluer Jab = 30.3, search. This is because at these redshifts the peak of the 
PPSNe spectra is shifted to ~ 40000 A, and observations at longer wavelengths 
represent an exponential increase in the observed flux. 

metal-free stars form at 1% of the observed total SFR density, which we model 
with a simple fit to the most recent measurements Q ^5.). In all cases we 
assume that A^|j^, the number PPSNe per of stars formed, is 10"'^. 



Comparing these models with ongoing surveys, we find that significant limits 
on PPSNe can already be placed out to moderate redshifts. Given the area 
and magnitude limit of the the Institute for Astronomy Deep Survey (js^, for 
example, one can constrain very massive star formation to < 1% of the total 
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SFR density out to a redshift of 2. 

Future missions should improve these hmits greatly. In Fig. 4, we show the 
PPSN constraints that would be obtained from three possible realizations of 
the planned space-based Joint Dark Energy Mission (JDEM). In the upper 
panel we adopt a limit of Jab = 27.7, corresponding to the wide-field survey 
planed for the Supernova Acceleration Cosmology Probe (SNAP) ^ realization 
of JDEM. With 2 months of data, the alternative Dark Energy Space Telescope 
(Destiny) ^ realization of JDEM would cover 7.5 deg^ of sky with spectroscopic 
observations down to a similar limit. In the central panel, we take a limit 
of Jab = 30.3 which would correspond to the deep-field SNAP survey. In 
the lower panel, we consider a redder survey with a limiting magnitude of 
Kab = 28.7, as appropriate for two months of observations from the Joint 
Efficient Dark-energy Investigation {J EDI) ^ realization of JDEM. 

Although planned with completely different goals in mind, these supernova 
searches will also serve as fantastic probes of PPSNe. Depending on the re- 
alization, JDEM may be able to place stringent limits on the rate density of 
even the faintest PPSNe out to z ^ 4:. For more luminous models, PPSNe 
will be constrained out to 2; ~ 6 and beyond. Thus, while a definitive mea- 
surement of metal-free stars remains allusive, narrow-band, spectroscopic, and 
SNe searches are fast closing in on these objects. The future of primordial star 
detection looks bright. 

I would like to thank UC Irvine and the organizing committee for hosting this 
enjoyable and informative workshop. 
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